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ABSTRACT 

In the hot, dense plasma of solar and stellar interiors, Coulomb potentials are 
screened, resulting in increased nuclear reaction rates. Although Salpeter's approx- 
imation for static screening is widely accepted and used in stellar modeling, the ques- 
tion of screening in nuclear reactions has been revisited. In particular the issue of 
dynamic effects has been raised by Shaviv and Shaviv who apply the techniques of 
molecular dynamics to the conditions in the Sun's core in order to numerically deter- 
mine the effect of screening. By directly calculating the motion of ions and electrons 
due to Coulomb interactions, the simulations are used to compute the effect of screen- 
ing without the mean-field assumption inherent in Salpeter's approximation. In this 
paper we reproduce their numerical analysis of the screening energy in the plasma 
of the solar core and conclude that the effects of dynamic screening are relevant and 
should be included when stellar nuclear reaction rates are computed. 

Subject headings: Equation of state - nuclear reactions, nucleosynthesis, abundances 
- plasmas - Sun:general 



1. INTRODUCTION 

Under the extreme temperatures and densities of the solar core, the plasma is fully ionized. 
The free electrons and ions interact via the Coulomb potential 
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U(r) = -. (1) 
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In this Coulomb system, nearby plasma is polarized by each ion. When two ions approach with 
the possibility of engaging in a nuclear reaction, each ion is surrounded by a screening cloud. 
Thus each ion is attracted by the electrons and repelled by the protons in its partner's cloud. The 
combined effect of the particles in the screening clouds on the potential energy of the pair of ions 
is referred to as screening. This screening effect reduces the standard Coulomb potential between 
approaching ions in a plasma to a screened potential which includes the contribution to the potential 
from the surrounding plasma. The reduced potential enables the ions to tunnel through the potential 
barrier more easily, thus enhancing fusion rates. 



Salpeter I (|l954|) derived an expression for the enhancement of nuclear reaction rates due to 



electron screening. By solving the Poisson-Boltzmann equation for electrons and ions in a plasma 
under the condition of weak screening (interaction << k B T), Salpeter arrived at an expression for 
the screenin g energy that is equival ent to that of the Debye-Hiickel theory of dilute solutions of 
electrolytes (IDebve & Huckelll 19231) . 

^screen -, (2) 

Ad 

where the Debye length, Ad, is the characteristic screening length of a plasma at temperature T 
with number density n, 

(3) 
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Although Salpet er's approximation for screening is widely accepted, several pap ers over the 
last few decades (e.g. IShaviv & Shaviv II 19961 : ICarraro et al. II 19881 : IWeiss et al. 11200 1|) have ques- 
tioned either the derivation itself or the validity of applying the approximation to hot, dense, 
Coulomb systems like the p lasma of the solar core. Various work deriving alternative formulae 
for electrostatic sc r eening dCarraro et al. 1 1 19 88; Qp her & Qpher I l2000i : IShaviv & Shaviv 1 1 19961 : 
Shavchenko 1 1 19991: Lavagno & Ouarati I I2OO0I : iTsytovich 1 12000) were refuted in subsequent pa- 
pers (see iBahcall et al. 1120021 for a summary of arguments in Salpeter's defense). However, the 
question of dynamic screening remains open. Dynamic screening arises because the protons in a 
plasma are much slower than the electrons. They are therefore not able to rearrange themselves 
as quickly around individual faster moving ions. Since nuclear reactions require energies several 
times the average thermal energy, the ions that are able to engage in nuclear reactions in the Sun 
are such faster moving ions, which therefore may not be accompanied by their full screening cloud. 
Salpeter uses the mean-field approach in which the many-body interactions are reduced to an av- 
erage interaction that simplifies calculations. This technique is quite useful in thermodynamical 
calculations that rely on the average behavior of the plasma. However, dynamic effects for the 
fast-moving, interacting ions lead to a screened potential that deviates from the average value. The 
nuclear reaction rates therefore differ from those computed with the mean-field approximation. 



Shaviv & Shaviv I (|1996L 1 19971120001 1200 1|) used the method of molecular dynamics to model 
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the motion of charges in a plasma under solar conditions in order to investigate dynamic screening. 
The advantage of the molecular-dynamics method is that it does not assume a mean field. Nor does 
it assume a long-time average potential for the scattering of any two charges, which is necessary in 
the statistical way to solve Poisson's equation to obtain the mean potential in a plasma. Shaviv and 
Shaviv attribute the differences between their simulations and Salpeter's theory to dynamic effects. 
Since their claims have been met with skepticism, we have conducted independent molecular- 
dynamics simulations to confirm the existence of dynamic effects. 



2. NUMERICAL METHODS 



For consistency, we developed molecular-dynamics simulations using the same basic nu- 
merical methods as Shaviv and Shaviv. To quell the skeptics, we evaluated the validity of 
each assumption and approximation before implementing it in our own simulations. Previ- 
ous publications iMao et al. 1120041 ; Klao llzoosl : iMussack etahl I2OO6L [20071 : IMussack 1 120071) 
describe this process and show that we did not find any errors in the methods used by Sha- 
viv and Shaviv. Once we validated their techniques, we used our independently developed 
simulations to examine dynamic effects in screening. 

Our numerical simulation consisted of a 3-dimensional box with 1000 particles (half protons 
and half electrons) interacting via the Coulomb potential. The temperature and density of the solar 
core (T = 1.6xl0 7 K, p = 1.6xl0 5 kg/m 3 ) provided the velocity distributions and inter-particle 
spacing. We applied periodic boundary conditions and the minimum-image convention. In order 
to deal with the long-range nature of the Coulomb potential, we implemented a cut-off radius. With 
this method, particles separated by a greater distance than the cut-off radius were not included in 
the potential sums as explained in section 12.11 Quantum effects were included through the use of 
effective potentials as desribed in section 12.21 Detailed analysis of the use of a cutof f radius and 



effective quantum potential as well as tests of our simulations can be found elsewhere (|Mao et al. 
20041 : IMao 1120081 : IMussack et aTkood. I2OO7I : IMussack lEoO^ . 



2.1. Dealing with Long-Range Forces 

The long-range nature of the Coulomb potential introduces additional challenges to a molecular- 
dynamics simulation. For many of the materials commonly modeled using molecular dynamics, 
short-range potentials such as the Lennard- Jones potential of a simple fluid 



u(r) = 4e 
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are appropriate. This potential decreases drastically with the separation of the two particles. There- 
fore it can be truncated quite reasonably at a distance of r cut = (2.5 - 3)<x. The Coulomb potential, 
on the other hand, decays slowly with particle separation. Particles at much greater distances must 
be included in order to accurately compute the force on each particle. Ideally, the contribution 
of each of the infinite image particles should be included in the potential sums and the potential 
should never be truncated. Unfortunately this would take infinite computing time and is therefore 
impractical. Fortunately, mor e efficient techniques have been developed to compute long-range 
interactions (see, for example. iFrenkel & Smit 1120021 : lAllen & Tildesley 1119871) . 



In their simulations, Shaviv and Shaviv took advantage of the screening of the plasma to 
shorten the effective range of the potential. They assume that interactions between particles that are 
farther apart than a few Debye lengths will be screened by the surrounding plasma. We have tested 
the truncation of the potential and det ermined that this assumption is appropriate for our simulation 
(fMussack et al. 1120071 : LMussack Il2007|) . This enables us to truncate the interactions beyond a few 
Debye lengths, reducing the infinte-range potential to a manageable long-range potential. 



2.2. Quantum Effects 



When particles in a two-component plasma are closer than the thermal deBroglie wavelength 

h 



^2nm a k B T 



(5) 



the effects of quantum diffraction and symmetry become significant. The thermal deBroglie wave- 
lengths for our system are shown in Table [TJ 

The average interparticle distance in our model is < r >= .203as. Since A p «; < r >, a 
classical treatment is generally sufficient for the protons. However, it will be common to have 
two electrons that are separated by a distance less than A e , requirin g quantum mech a nical treat- 
ment. We therefore include the quantum approximations used by 



the effective pair potentials derived for a hydrogen plasma by iBarker 



Shaviv & Shaviv 



(1996b, i.e. 



dl97lh : lDeutschNl977h : 





Mass 


A« 


Particle 


(m e ) 


(«b) 


Proton 


1836 


.008488 


Electron 


1 


.3637 



Table 1: Thermal deBroglie Wavelenghts at T = 1.6jc10 7 K 
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Deutsch et al. I (|l978lll979|) . Quantum diffraction effects are described by 



(d) _ e » e f3 

V/? - _ 



where 



and jj. a a is the reduced mass of the pair. 



1 - exp ( — 
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The exclusion principle between electrons is included by adding 



v ( £ = k B T In 2 exp 



A In 2 



(9) 



These potentials only differ from the classical Coulomb potential at short distances, where quan- 
tum effects are relevant. Fig. Q] shows the contribution to the effective potential from quantum 
corrections to the electrons and protons. It is clear that quantum effects for electrons can not be 
neglected at the high temperature of the solar core. The correction for proton-proton pairs is neg- 
ligible, but electron-proton pairs do experience diffraction effects at close distances. Therefore we 
include quantum corrections for both electron-electron and electron-proton interactions. 



2.3. Method 

The goal of this work is to compute the energy exchanged between a pair of approaching 
protons and the surrounding plasma in order to evaluate the screening energy and the dependence 
of the screening energy on the relative kinetic energy of the pair. In order to accomplish this goal, 
we use a method similar to that of Shaviv and Shaviv. The technique for tracking approaching 
proton pairs and the energy they exchange with the plasma is described below. 

1. For each proton i, find all protons approaching it. Choose the closest approaching proton 
and call it the partner of proton L 

2. Track the pair as they approach and then move apart. Record their distance of closest 
approach, r c , and the energy of the pair at their closest point 

^ air = kinetic + 7- (10) 



-6- 



1.50 




0.00 1 1 1 1 1 1 

0.05 0.10 0.15 0.20 0.25 0.30 



Fig. 1 . — Comparison between effective pair potentials including quantum effects and the classical 
Coulomb potential (solid line) in a plasma under the conditions of the solar core, ^(triangles) 
and V^squares) include both diffraction effects and symmetry effects, while V ep (circles) includes 
only diffraction effects. 
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3. Continue tracking the pair until they are separated by a distance Rf. For our simulations, 
Rf = 2.0a B is a sufficient distance to represent the separation of the pair. At this point, record the 
far-apart energy of the pair 

E f = E f + — (11) 

pair kinetic Rf' 

Beyond this point, the pair is no longer tracked. 

4. Find a new partner for proton i, and repeat the process for i and its new partner. 

Since a proton pair's departure from their closest point to a distance Rf is symmetric to an 
approach from Rf to r c , we use the information obtained from the departure to examine the dynamic 
effect of the plasma on the screening enhancement for interacting proton pairs as they approach a 
separation at which nuclear reactions are possible. The dynamic information of each proton and 
its partner du jour are recorded at each time step. 



3. RESULTS 

First we examine the distribution of the distance of closest approach of proton pairs. Fig. [2] 
shows the number of pairs that achieved each closest distance r c . We see that most pairs are able 
to get at least as close as the average distance of neighboring protons < r p >= 0.256cib- As seen in 
Fig. [31 pairs of protons with higher far-apart kinetic energy E^-^ are more likely to get closer to 
each other than pairs with lower far-apart kinetic energy. 

The screening energy is the energy exchanged between a pair of protons and the surrounding 
plasma. It represents the combined effect of the neighboring protons and electrons on the inter- 
acting pair. Salpeter's mean-field theory gives an expression for the average screened Coulomb 
potential that includes the energy exchanged between interacting protons and the plasma during 
their approach 

l(r) = -{\-e- rlx °) . (12) 



4 screen V 

r 



Salpeter's theory averages the effect of the kinetic energy of each particle so that the screening 
energy only depends on distance. However, our simulations show that the screening energy also 
depends on the kinetic energy of the protons. Protons with high kinetic energy tend to gain less 
energy from the plasma during the interac t ion, while p rotons with low kinetic energy gain more 



energy from the plasma. IShaviv & Shaviv I (|2000l . l200l|) refer to this as the dynamic effect. 



For the purpose of comparison, we define the screening energy in the simulations to be the 
difference in the energy of the pair when they are far apart and when they are close together. 



(13) 
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e 1 

^pair = ^kinetic + ~ (^) 
' r 



? 2 

'pair ^kinetic ' 

Since all pairs have the same far-apart distance Rf, 

? 2 

' potential 



nair kinetir n 

(15) 



EjLnttal = = C ° nStant ( 16 ) 

This yields a screening energy of 

^screen - ^kinetic _ ^kinetic + ~ ^potential 1 (^) 

' c 

A key ingredient in this equation for the screening energy is the difference in the kinetic energy 
of a pair when the partner protons are far apart and when they are at their closest separation. 

A£kinetic = (^kinetic _ ^kinetic) 0^) 

Here we examine the relationship between the average change in the kinetic energy of a pair and 
their distance of closest approach, shown in Fig. HI This energy change is compared with the 
Coulomb potential and the screened Coulomb potential from Debye-Hiickel theory. We see that at 
the distance of closest approach, the average kinetic energy exchanged between a pair of protons 
and the surrounding plasma has the form of the screened Coulomb potential. 

<AE kinetlc >l c ~ f-g-'cMi. (19) 
r c 

On average, the protons interact in a mean field potential of 

e 1 

V(r) = —e- rlAD (20) 
r 

as Salpeter's theory states. 

We can rewrite the average screening energy as 

„2 



^screen >lr c ~ (l 6 l/ d J ^potential- @1) 



As R f -> oo, £^ otential -> 0, leaving us with 



*^ ^screen ^ I r c — ^screen 0"c)' (22) 
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Fig. 4. — Dependence of the average kinetic energy change <A£ , k i ne tic> on the closest distance r c , 
compared with the Coulomb potential and screened Coulomb potential. 
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The average screening energy in the simulation is the same as the average screening energy derived 
in Salpeter's theory. This confirms that on average, the system can be described by the screened 
Coulomb potential of the Debye-Hiickel theory as Salpeter contends. 

In order to see the dynamic effect on screening, we must look at the relationship between the 
screening energy and the far-apart kinetic energy of proton pairs. Fig. |5] shows the average energy 
gained from the plasma by pairs of protons with a given far-apart kinetic energy in the simulation. 
For comparison, the Debye-Hiickel screening energy computed at a distance equal to the average 
closest-approach disance of pairs of protons with each far-apart kinetic energy is also shown. The 
screening energy of a pair of interacting protons depends on the relative kinetic energy of the pair. 
Pairs with far-apart kinetic energy greater than the thermal energy tend to gain less energy from 
the plasma than the mean-field screening energy, while pairs with lower far-apart kinetic energy 
gain energy more energy from the plasma than the mean-field screening energy. 



^screen -^screen fof E^^^ < kjjT (23) 
^screen ^ ^screen fOT ^jjjjjgfjQ - > kgT (24) 

When the screening energy from the simulation is averaged over pairs with all E( inetic , the 
result is the same as in Salpeter's theory. This reaffirms that Salpeter's theory is appropriate as a 
mean-field treatment. 

However, the dynamic effect becomes relevant when pairs with high relative kinetic energy 
are singled out, as in solar nuclear reactions. Pairs with far-apart kinetic energy near the Gamov 
energy of p+p reactions (4.8k B T) will, on average, lose energy to the plasma rather than gaining 
energy. Because of this dependence of the screening energy on the relative kinetic energy of a pair 
of interacting particles, we conclude that the mean-field approach is not appropriate for computing 
nuclear reaction rates. 



4. CONCLUSIONS 

Our simulations confirm that the static Debye screening potential describes the average be- 
havior of the plasma as a whole. Therefore it remains the appropriate treatment for mean-field 
properties of the system. However, dynamic effects are important when particles are singled out 
based on their velocities. Pairs of particles with greater relative velocities experience less screening 
than pairs of particles with lower relative velocities. This effect is relevant in nuclear reactions in 
the Sun because they occur at energies that are several times greater than the thermal energy. The 
fast protons from the high-end tail of the Maxwell-Boltzmann velocity distribution are most likely 



-13- 




Fig. 5. — Dependence of the screening energy from the simulation ^screen (squares) on the far- 
apart kinetic energy e{. . . The Debye-Hiickel screening energy computed at the averege closest- 
approach distance of pairs of protons with a given far-apart kinetic energy is shown (circles) for 
comparison. 
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to take part in these reactions. Thus the mean-field approach is not valid in this case. We agree 
with Shaviv and Shaviv that the formalism for screening in solar nuclear reaction rates should be 
amended to account for this dynamic effect. 

We thank Rajiv Kalia, Aiichiro Nakano, and Priya Vashishta for their guidance on the molecular- 
dynamics simulations. This work was supported by grant AST-0708568 of the National Science 
Foundation. 
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